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Abstract 

Issues of proton decay are examined in supersymmetric standard model 
with U(l) flavor symmetry. Dimension five proton-decay operators which 
arise generically are controlled by the flavor symmetry. We show that unlike 
the minimal supersymmetric SU(5) case the proton decay modes containing 
charged lepton can have large branching ratios if the dimension five operators 
of left-handed type dominate. Measuring the branching ratio of the electron 
mode to the muon mode may reveal the mechanism of the neutrino mass 
generation. The case with vanishing charges for Higgs doublets marginally 
survives the present experimental bound on proton life time. 
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I. INTRODUCTION 



With particle contents and gauge symmetry of the standard model, nucleon decay may 
take place through non-renormalizable interaction. In fact, the gauge invariance allows 
dimension six operators with baryon and lepton number violation. They may be induced 
by exchange of X, Y gauge bosons and colored Higgs multiplets in a grand unified theory 
(GUT), or may just arise as non-renormalizable operators at a fundamental scale which 
we assume to be near the Planck scale. The operators of the latter origin are, however, 
suppressed by the Planck mass squared, and thus the nucleon decay induced by them is 
suppressed to a level where near future experiments will not be able to reach. 

The situation is very different if one considers the supersymmetric (SUSY) extension 
of the Standard Model. Existence of superparticles allows gauge invariant dimension five 
operators which can induce nucleon decay after superparticle dressing [[]]].[] These operators 
are very dangerous because they are suppressed only by a single power of the Planck mass. In 
fact, for the superparticle masses around 1 TeV, present proton decay experiments constrain 
the mass scale of the dimension five operators much larger than the Planck mass, or in other 
words, their coefficients should be much smaller than unity when normalized by the Planck 
scale. This is indeed embarrassing if one believes the widely accepted argument on the 
generality that all operators which are allowed by symmetry should arise with order one 
coefficients. There are many attempts to explain the smallness of these generic dimension 5 
operators. They include 1) imposing some symmetry such as the family symmetry JJ,[5],[S]|)|] , 
the discrete gauge symmetry 0, the Peccei-Quinn symmetry @,|], the U(l)yi symmetry [TIJ 



and the R symmetry |TT],|T2|, and 2) attributing to configurations of quarks and leptons in 
extra spatial dimensions |T3|JT^ . 



In this paper, we shall re-examine the case of U(l) flavor symmetry JT5[ and closely 
investigate how the nucleon decays can be (or cannot be) suppressed in the context of 
the minimal supersymmetric standard model (MSSM). Furthermore, we will explore the 
correlation between the branching fractions and the family structure [0-ffl. 

We will show that, unlike a typical SUSY SU(5) GUT, the proton decay modes containing 
charged lepton can have sizable branching ratios and the ratio of the [i mode and the e mode 
may give us a crucial information on the neutrino mass generation. We will also study the 
proton life time. Keeping in mind that our argument based on the flavor symmetry contains 
uncertainties in coefficients of operators, we will argue that the case with vanishing charges 
for the Higgs doublets marginally survives the present proton decay bound. The constraint 
is relaxed for negative Higgs charges. We will briefly discuss how to generate higgsino mass 
/i in this case. 



We assume .R-parity conservation, and thus we do not consider the most dangerous dimension 
four operators. 
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II. U(l) FLAVOR SYMMETRY AND CHARGE ASSIGNMENT 



A hypothetical flavor symmetry is a symmetry to explain the Yukawa structure. Here 
we consider the Froggatt-Nielsen mechanism [15[] in which we assign Froggatt-Nielsen U(l) 
charges q to the MSSM superfields Q and —1 to the Froggatt-Nielsen field X. The U(l) 
symmetry is spontaneously broken when the X field develops a vacuum expectation value. 
We define A as A = (X)/M p \, with M p \ being the Planck scale. In the MSSM, we have the 
following superpotential 



W = y l SQiUfH u + y%QiDfH d + y % £EfLjH d 



(1) 



where i,j = 1,2,3 are generation indices, etc. are Yukawa couplings, whose magnitudes 
are governed by the U(l) symmetry 



Vv 



-u a .+h u 



Vd 



ij \ qi+dj+h d 



y% = ./ /:• A' ' ' ' 



(2) 



Here f$ etc. are somewhat arbitrary constants which are typically of order unity. In order 
to obtain physical masses and mixings, we must translate fields from a flavor basis into a 
mass basis through unitary matrices as 



UfvuUu 



diag(y U} y c ,y t ) 
U^yoU D = diag(y d ,y s ,y b ) 
UeVlUl = diag(y e ,y^y T 



(3) 
(4) 
(5) 



from which we obtain the CKM matrix as Vkm = UqUq. 

Following the conventional wisdom, we identify A with the Wolfenstein parameter A ~ 
0.22, and determine the U(l) charges of the MSSM fields. The experimental values of the 
masses and mixings near the Planck scale are approximated as 



V,. 



A, V cb ~ A 



V, 



ub 



(6) 



and 





: m c 


■ m t 


~ A 8 


A 4 


1 


m d 


■ m s 


m b 


~ A 4 


A 2 


1 


m e : 




m T 


~ A 5 


A 2 
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In this paper, we will consider the following class of charge assignments: 



9i 
h 



q 2 = q 3 + 2, 
u c 2 = u c 3 + 2, 



93 + 3, 
e c 3 + r- 

d% + 1, d c 2 = dl 
l 3 + n + m, l 2 = h + n, 



+ 2 



r 

l z ■ 



(7) 
(8) 
(9) 



(10) 

(11) 
(12) 

(13) 
(14) 



where n, m,p and r are integers with < n < 2,0 < m < 5 — n,3 < p < 5 and 4 < r < 
5. Here, eq. (|i~0D is dictated by the magnitudes of the CKM matrix elements, eq. (|i~l~D and 
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eq.(|13D are obtained from the up- and down- type quark mass ratios and the mixings, and 
eq. (|i~2"D and eq.(14) are determined by the charged lepton mass hierarchy. Uncertainty of 



the coefficients f}j etc. allows some ambiguities of the charges, which are parameterized by 



p and r. Furthermore it follows from the fact that 



m. 



and y t 



+ l 3 , q 3 + u c 3 + h u = 0. 



(15) 



The parameters m and n are determined by the oscillation of the atmospheric |T7] and 
solar [T3|] neutrinos. Here we assume that the see-saw mechanism |B| operates to generate 
neutrino masses and mixings. The overall scale of the neutrino masses is fixed by the scale 
of lepton number violation, and the ratio of the masses as well as the mixing angles are 
controlled by the U(l) charges. We will consider the following two cases, both of which 
reproduce the bimaximal mixing angles: 
i) the lopsided type [ED 



n = 0, 



m 



(16) 



and 



ii) the anarchical type [|21|j22 



n 



0. 



m 



0. 



(17) 



with the right-handed neutrinos having appropriate U(l) charges. 

To make our argument explicit, we will consider the following two types of charge as- 
signment. 

1) best-fit charge assignment 

In the first type of charge assignment, we take 



V 



■5. 



n 



0. 



or 1, 



(18) 



in eqs. (pT0|) — (|14D . We call this assignment the best-fit one. Notice that it is not consistent 
with the GUT symmetry where the Q and U , for instance, are in a single multiplet and 
thus have a common U(l) charge. 
2) GUT -inspired charge assignment 

The second type of charge assignment we will discuss is the one consistent with the SU(5) 
GUT. Specifically we consider 



d\ 



qi = Q3 + 3 

di + 1, 



q-2 
d c 3 



<?3 + 2, 

(n = 0, 



(P = 3, r = 4), 



m 



(19) 
(20) 



The charge assignment for u c would predict the up-quark Yukawa coupling of A 6 , much larger 
than the actual value. Thus in this scheme, the small up-quark mass should be attributed to 
either accidentally small coefficients fy of order A 2 or the cancellation between the matrix 
elements. 
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III. PROTON DECAY 



Now we would like to discuss proton decay. Gauge invariance allows the dimension five 
operators in the superpotential, which break the baryon number as well as the lepton number 
conservation and are divided into an LLLL part and an RRRR part: 

W = ^CgVQiQjQkLt + ^Cf l EfUfU°Df. (21) 

Here CT and Cg w are some coefficients, and M denotes the fundamental scale of the 
theory, which we identify with the Planck scale M p \. If these coefficient were of order unity, 
they would cause proton decay with life time several magnitudes shorter than the present 
experimental bound. In the framework of the Froggatt-Nielsen U(l) symmetry, we expect 
suppression 

C ijkl ^ ^qi+qj+qk+h^ yg"y*'A~^~ hd A*~ tt 5 +,J ~ < ?), (22) 

C% kl ~ \ e i+ u 1 +u % +d i(~ y%y k l\- h ™- h *\ e i- qi + u t-' l <°). (23) 

Here we have omitted constants in front of the powers of A. It is interesting to note here that 
the proton decay is suppressed if the sum of the two Higgs charges is negative, h u + hd < 0. 
We will come back to this point later on. 

After wino and higgsino dressing,^ we obtain proton decay amplitudes |24H29| f\ 



Amp(p - K + u e ) ~ -L(A\ + AH (24) 
Ampip^K+^^L^l + A^ (25) 

Amp(p -> K+v T ) ~ -L(Al +Al+A 3 + A 4 ), (26) 



where 



2 We assume that flavor mixings in the quark-squark-gluino interactions are small enough to avoid 
the dangerous flavor changing processes, and simply ignore gluino dressing diagrams in this paper. 
If these generation mixings are sizable, the gluino dressing diagrams give considerable contributions 



to proton decay [23|. 



3 Proton decay to pion suffers from Cabbibo suppression, and thus we will not consider it in this 
paper. 

4 Cf kl in eqs. fl27|) and @ should be replaced by C™ npq (U%) m \U%) n j (U%) p k (U L ) q l . However, 



since the off-diagonal elements of the mixing matrices are suppressed by some powers of A (includ- 
ing A ), we deduce that C?£ nm (U%)J(U§) n j (U%) p k (U L ) q l ~ C^ 1 in this framework. The same 
argument is applied to Cr. 
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A\ ~ g 2 2 V cd V cs Cf l \ (27) 

^ ~ glv td v ts cf lfc , (28) 

Ag-M/.V^Cf 11 , (29) 
A, ~ y t y T V u C™ 12 . (30) 

In eqs. ( 2~T|) and (E§), we have omitted other contributions proportional to Cj} 2k , C}} 3k , • • ■, 



which are of the same order as the terms explicitly written there. In eqs. - (|26|) / 

represents a loop factor, / ~ l/(167r 2 msusY), where m SU sY is a representative sparticle 
mass scale. Degeneracy of the squark and slepton masses in different generations is implic- 
itly assumed in the above evaluation. That is, we do not consider the so-called effective 
supersymmetry where the squarks and sleptons in the first two generations are very heavy. 

The amplitudes of the decay n — > K°u are the same as those of p — > K + u up to factors 
of order unity. 

Using a numerical relation 



VtVr = 2m2 ' ^ gl (^tan/5+ ~ -A 2 fl£tan/3, tan/3 > 3 (31) 



at the Z-boson mass scale, we can rewrite the coefficients A as 

A\ ~ A\ ~ gly t y b \- h ^- hd+s+s+ ^- h \ (32) 
As ~ A A ~ g 2 2 yty b X~ h - hd - S+5+pt -^-, (33) 



where 6 = 2q 3 — u\ — e\. Notice that the S dependence of the LLLL operators and that of 
the RRRR operators are inverse. 

As for charged lepton modes, we obtain 

Amp(p - K°e + ) ~ La\, (34) 

Amp(p -> #V) ~ ^ (35) 

where 

4 = <fe 2 E Ci 1 ^* ~ gly t y b \~ hu ~ hd+5+&Jt '^ k ~ 1 ^ . (36) 

i 

Since the proton cannot decay to the tau lepton, the decay amplitudes do not contain the 
third generation down- type Yukawa coupling constant. Therefore, contributions from the 
RRRR operators to the charged lepton modes are negligible. 



A. Branching Ratios 

Let us first argue the branching ratios of the proton decay. When the LLLL operators 
dominate in the decay modes to the neutrinos, we find from the above consideration 
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Thus, we predict that 



r(p -> K°e+) (n+m) r(p -> gV) 2n 

r(p — " ' r(p " ' 1 J 

where T(p — > = J2k r(p — > K + Vk) is the decay width into three types of the neutrinos. 

In particular, for n = 0, which is strongly suggested by the large mixing between v T and 
we conclude that the decay rate of the muon mode is comparable to that of the neutrino 
mode. Furthermore, the decay rate to the electron is controlled by the U(l) charge of the 
li, and thus we may be able to distinguish the anarchical type assignment for the leptons 
from the lopsided one by measuring the ratio T(p — > K°e + )/T(p — > K Q ^i + ). 

This consequence is a striking contrast to the case of the minimal SU(5) SUSY GUT in 
which the proton decay operators are induced by exchanging the colored Higgs multiplets. 
In fact, in the SU(5) case the charged lepton modes are negligibly suppressed |25| , |26|j . This is 



because the decay amplitudes into the charged leptons are proportional to the very small up- 
quark mass, while those into the neutrinos can be proportional to the charm- or top-quark 
mass. 

On the other hand, when the RRRR operators are dominant, the proton mainly decays 
through p — ► K + v T . Thus we expect 

Which of the operators, LLLL or RRRR, dominates depends on the charges one assumes. 
When we take the best- fit charge assignment, I2 — h = and p — 5, whereas 5 is a parameter 
which is not fixed. Then, for 

5< 0.9-0.31n(tan/3/3), (40) 

the LLLL contributions dominate over the RRRR ones. 

In the GUT-inspired case, I2 — h = 0,p = 3 and 5 = 0.Q It follows from eq.(^) that for 
small tan f3, the LLLL amplitudes are of the same order as the RRRR amplitudes. Thus, 
we expect that the muon mode is comparable to the neutrino mode and the electron mode 
is suppressed by A 2 . Namely, 

r <r*"f)~l, r ^ A '° e+ >~A* (tan/3. 3). (41) 

On the contrary, for large tan f3, the RRRR contributions are enhanced and thus the neutrino 
mode, p — > K + u, dominates. 



5 If the theory is really embedded into a GUT group, one has to take into account the contribution 
from the colored Higgs exchange. In the GUT framework, our argument given here is valid when 
some mechanism operates to suppress the proton decays mediated by the colored Higgses and thus 
the genuine dimension five operators controlled by the flavor symmetry are the dominant sources 
of the proton decay. The suppression mechanisms of the proton decay operators coming from the 
exchange of the colored Higgses have been proposed, for instance, in Refs. |^,10|,30,31|. See also 
Ref. [32 1 and references therein. 
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B. Decay Rates 



Next we will derive constraints on the charge assignment of the MSSM fields from the 
null results of the proton decay searches. In the framework we are considering, the decay 
rates for the modes p — > K + u, K°e + , K°p, + and n — > K°v can be comparable in magnitudes. 
However, since the most severe experimental constraint comes from the decay mode p — > 
K + u, all we have to do is to study this mode. 

The partial decay rate for the mode p — > K + v is calculated as 

F(p - K + u) ~ ^ ( 1 - -1 |i? L /3 p Am m + ^^Ampfl^rl , (42) 



32tt V ml J II 



k 



where 



Am PLk ~ ^^A-^-^ +5+8+fe -' 3) , (43) 
Ampi, ~ ^^ 6 A-^-^- 5+5+ ^. (44) 

Here, m p and are the masses of proton and K + respectively, f w is the pion decay 

constant, Rl,r represent renormalization effects of LLLL and RRRR operators from M to 
1 GeV, and a p and f} p are the hadronic matrix element parameters. Hereafter, for simplicity 
we fix the parameters as msusY = 1 TeV, \a p \ = 0.015 GeV 3 , |/3 P | = 0.014 GeV 3 [^3f,i?L = 
10.2 and R R = 6.5. f\ 

For the best-fit charge assignment, in order not to conflict with the experimental bound 
r(p — > K + v) > 1.9 x 10 33 yr [34], the U(l) charges must satisfy the following inequalities: 

yty b \- hu ~ hd+5+8 £ 4.2 x io- 8 , (45) 
y t y b \- K - hd - 5+wt -^ < 6.2 x io- 8 . (46) 

Substituting the approximate values m t ~ 110 GeV and m;, ~ 1.0 GeV near the Planck 
scale, the above are rewritten as 

\- hu - hd+s tan f3< 2.1, (47) 

x - hu - hd - s tetf£ < 64 _ (4g) 



Thus, we obtain 



K + h d < -0.2 + 5-0.7 ln(tan/3/3), (49) 
/in + ^^l-8-5-1.31n(tan/3/3). (50) 



Rl,r are evaluated by solving one-loop renormalization group equations due to gauge interactions 
from a high energy scale (~ 2 x 10 16 GeV) to 1 GeV. Inclusion of the Yukawa couplings and effects 
above this high energy scale do not qualitatively change our results. 
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The proton decay constraint is relaxed most when tan/3 ~ 3 and 5 = 1, which gives the 
upper bound 



h u + h d <l, (51) 

and hence we find that the case h u + h d = survives the proton decay constraint for very 
low tan/?. On the other hand, for larger tan/3, the sum of the two Higgs charges must be 
negative for any choice of 5 to satisfy the proton decay bound. 

For the GUT-inspired case where 5 = and p = 3, the constraints on the Higgs charges 
are modified as 

h u + h d < -0.2-0.71n(tan^/3), (52) 
K + h d < -0.3-1.31n(tan^/3). (53) 

Enhancement of the RRRR contributions arises from the fact that difference between u\ and 
«3 is smaller than that in the best-jit case. Thus we conclude that the charge assignment 
h u + h d = is marginally allowed for very small tan/3. 

In deriving the constraints above, one has to keep in mind that there are some uncer- 
tainties in evaluating the proton decay rate. They include 

1. The coefficients of the proton-decay dimension five operators which have been set order 
unity may be accidentally small, like m u in the GUT-inspired case. 

2. There are some uncertainties in a p and f3 p . Taking the smallest allowed values for 
them, the proton decay amplitudes become five times smaller than our estimate. 

3. We have taken the representative sparticle mass rrisusY — 1 TeV. However, it becomes 
larger and the proton decay gets suppressed if we adopt smaller wino and higgsino 
masses and larger squark and slepton masses. 

These could reduce the proton decay rate, and then the constraints for the Higgs charges 
obtained above would become somewhat relaxed. 



On the other hand, negative charge assignment for the Higgs doublets |30[ survives the 
proton decay constraint in a wider region of the parameter space as the proton decay rate 
is reduced by \- 2 ( h ^+ h d)_ 



IV. fx PARAMETER 

Here we would like to make a brief comment on the higgsino mixing parameter /i when 
h u + h d is negative. In this case, neither the term [iiH u Hd]e 2 nor \X j 'M p \) n H U H d ] g 2 (with 
n > 0) is invariant under the Froggatt-Nielsen U(l) symmetry. Thus the generation of 
the n term is somewhat contrived. Consider the Giudice-Masiero mechanism p5] with 
the term [ZyM pl (XyM pl )-^ +h ^H u H d ] 8 2 q2. This gives the fi term of F|/M pl A- ( ^ +hd) ~ 
ms/2^ hu+hd ^- This is much smaller than the electroweak scale unless the gravitino mass 7713/2 
is large. This is, in fact, the case in, e.g. anomaly mediation where 1713/2 is around a few 
tens TeV. Thus as far as \\ hu+hd \ > 10~ 2 , we obtain the weak scale fj, parameter. However, 
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the Higgs mixing parameter, B, is not suppressed and in general as heavy as m 3 / 2 . Then 
we will have difficulty in obtaining the correct electroweak symmetry breaking. 

A possible mechanism is to introduce a singlet S with a positive U(l) charge so that 
[SH u H<i]q2 is invariant under the flavor symmetry. It is assumed that S field does not have 
a VEV in the SUSY limit. SUSY breaking generally generates a VEV of S, which is nothing 
but the fi | 36|j . An explicit realization was given in Ref. | 37| , where the charges of relevant 
fields should be adjusted to obtain the \x parameter of correct order of magnitude. Thus the 
choice of a negative h u + hd does not immediately cause trouble to generate the n parameter, 
though the way to do is rather restricted. 



V. CONCLUSIONS AND DISCUSSION 

In this paper, we have considered the proton decay when the dimension five operators 
are controlled by the U(l) flavor symmetry. We have shown that the charged lepton decay 
modes can have sizable branching ratios when the LLLL operators dominate over the RRRR 
ones. In this case, measurement of the branching ratio of the e mode and the \i mode will 
provide an important information on the charges for the doublet leptons and thus reveal 
how the neutrino masses are generated. On the other hand, if the RRRR operators are the 
main sources of the proton decay, the neutrino mode dominates in the proton decay. We 
have also investigated the proton decay rate and showed that the proton-decay suppression 
due to the U(l) flavor symmetry marginally survives the present proton decay bound even 
when the Higgs doublets do not have the U(l) charges. The bound is less severe, if we take 
negative charges of the Higgs doublets. 

We have made several assumptions in deriving the conclusions drawn above. Firstly 
we have assumed that non-renormalizable operators are suppressed by the four- dimensional 
Planck scale. However, there are many models in which the fundamental scale is lower than 
the Planck scale. In this case, it is natural to expect that the proton-decay dimension five 
operators are suppressed only by this low fundamental scale and hence the proton decay rate 
is considerably enhanced. To survive the present experimental bound, one needs to assign 
the negative U(l) charges to the Higgs doublets in this case. Arguments on the branching 
ratios remain unchanged. 

Another important assumption we have made is the degeneracy of the squark and slepton 
masses in different generations. Let us consider the case where the superparticles of the first 
two generations are heavy, which is realized when an anomalous U(l) flavor symmetry also 
mediates SUSY breaking. In this scheme, as far as the scalar tops I^r and the right-handed 
scalar tau ?r are light like in the case of GUTs, the amplitudes of the LLLL operators are 
greatly suppressed while those of the RRRR operators remain unsuppressed |38|.|39| . Thus, 



one finds that p — > K + v tends to dominate over K°e + , K°fi + modes. On the contrary, when 
only Il^r are light, the RRRR operators are suppressed to the same extent of the LLLL 
operators. Therefore, the discussions about the branching ratios are not altered. 
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